Perturbations on steady spherical accretion in Schwarzschild geometry 
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The stationary background flow in the spherically symmetric infall of a compressible fluid, coupled to the 
space-time defined by the static Schwarzschild metric, has been subjected to linearized perturbations. The 
perturbative procedure is based on the continuity condition and it shows that the coupling of the flow with 
the geometry of space-time brings about greater stability for the flow, to the extent that the amplitude of the 
perturbation, treated as a standing wave, decays in time, as opposed to the amplitude remaining constant in 
the Newtonian limit. In qualitative terms this situation simulates the effect of a dissipative mechanism in the 
classical Bondi accretion flow, defined in the Newtonian construct of space and time. As a result of this approach 
it becomes impossible to define an acoustic metric for a conserved spherically symmetric flow, described within 
the framework of Schwarzschild geometry. In keeping with this view, the perturbation, considered separately as 
a high-frequency travelling wave, also has its amplitude reduced. 
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I. INTRODUCTION 

The spherically symmetric model of astrophysical accretion 
continues to enjoy an abiding appeal among researchers in ac- 
cretion astrophysics, starting with the seminal paper written 
by Bondi 1 1] more than half a century ago. The basic simplic- 
ity of this model notwithstanding, it is actually quite appro- 
priate for many realistic aspects of accretion processes, and 
certainly more than anything else, this model also allows a 
clear insight to be had into the related physics, much of which 
is often of quite an involved nature. As a result the spherically 
symmetric model is frequently the starting point from where 
it becomes possible to devise theoretical models of increasing 
physical complexity for accretion processes. 

While general questions related to astrophysical accretion 
have been addressed and studied from various perspectives all 
along, it was not too long before Bondi's original treatment, 
carried out within the Newtonian construct of space and time, 
was extended to embrace a general relativistic description of 
spherically symmetric accretion. An early work in this regard 



'Electronic address: tptn@mahendra.iacs. res. in' 
^Electronic address: tpnc@mahendra.i acs.res.in| 
tElectronic address: tpjkb@mahendr a.iacs .res. in] 
§Also at: Positional Astronomy Centre 
P-546, Block 'N' New Alipure, Kolkata 700053, India 
^Electronic address: ,iakr@iucaa.emet.in.i 



was reported by Michel f^, which was followed by a spate 
of later works, of which some took up various issues ranging 
from the stability of solutions to their exact nature and critical 
aspects iliiilTlii. 

Stability of spherically symmetric accretion has been stud- 
ied long and studied variously il[l^[IIl[ll[ll[ll[ll. With 
especial regard to the stability of spherical accretion on to a 
black hole, Moncrief's study |6] has shown no evidence of 
the development of any instability under the influence of a 
linearized perturbation on the standing background flow. The 
treatment presented in this paper returns to the same general 
theme, but the methods applied here, and the motivation that 
has prompted them, are different. First of all, as opposed to 
Moncrief's approach of perturbing a scalar potential, whose 
gradient is prescribed to be the velocity of the ideal fluid, in 
this paper the perturbation scheme that has been adopted is 
centred around the continuity condition. This follows the pro- 
cedure employed earlier by Petterson et al. fTl ] and Theuns & 
David 1 12], respectively, in their studies of spherical accretion 
in the non-relativistic regime. The stationary solution of the 
continuity equation gives a first integral, which, within a con- 
stant geometric factor, is actually the matter inflow rate. The 
full description of the flow will imply that complete solutions 
of two coupled fields — the velocity field and the density field 
— will have to be obtained. Both these fields are connected 
to the matter accretion rate, and, therefore, perturbing the ac- 
cretion rate about its constant stationary value will lead to a 
wave equation for a single perturbed field that will also convey 
enough useful information on the stability of both the velocity 
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and the density fields. 

While this is the general procedure that has been followed, 
the primary objective of the whole exercise has been to see 
whether or not through this perturbative technique it should 
be possible to establish an acoustic geometry for a flow that 
is at the same time described under fully general relativistic 
conditions. In recent years fluid analogue gravity has been a 
subject that has, extending over diverse types of fluid flows, 
called much attention upon itself 11151 

IllllllllllllMlllllilllllll UM- In fluid dynamical pro- 
cesses a critical point is the point where the speed of the bulk 
motion matches the speed of information propagation through 
the fluid. Depending on the direction of the flow, information 
of any event occurring in either the super-critical region of the 
fluid (where the bulk flow is faster than the speed with which 
any information can travel) or its sub-critical region (where 
information propagation overrides the bulk flow) will not per- 
colate into the other region, through a surface determined by 
the critical condition. For an ideal fluid, therefore, the critical 
condition defines a barrier, which can be viewed as the fluid 
analogue of the event horizon of either a black hole or a white 
hole, according to the direction in which the flow proceeds. 
It can be easily appreciated that for spherically symmetric ac- 
cretion solutions the sonic surface defines the event horizon 
of an acoustic black hole. This feature has been studied ex- 
tensively and understood well by now in both the Newtonian 
framework | [T5|, IztI Isoll as well as in the general relativistic 
framework 11221 um . The latter context is very interesting be- 
cause it combines the metric properties of curved space-time 
with similar properties of the fluid itself flowing in the same 
space-time. 

This work purports to investigate the same analogue geo- 
metric properties of the flow, described in the Schwarzschild 
metric. A study along these very lines has been reported ear- 
lier by Bilic |20], whose work, however, is based on the usual 
practice of studying a perturbation on a scalar potential func- 
tion. In contrast, as it must be emphasized once again, the 
present work approaches the whole question from the view- 
point of the continuity condition in the flow. This difference 
of approach turns out to be significant because the result deriv- 
ing from the latter line of attack is negative, quite unlike what 
has been shown through the method used by Bilic [20] . This is 
rather curious because in the non-relativistic framework both 
paths can be shown to lead to the same end [ TsilZTilsofl . It has 
been seen here that the coupling of the flow with the spheri- 
cally symmetric geometry of space-time acts like a dissipative 
effect, and breaks down the Lorentz invariance that should be 
necessary to define an analogue metric. The invariance is re- 
stored in the Newtonian limit. By way of comparison and 
by following the same mathematical prescription, one could 
invoke a similar feature that arises because of viscosity in 
shallow-layer incompressible fluid flows where information 
propagates as gravity waves [21]. When viscosity is made 
to vanish, it becomes possible to establish an analogue black 
hole (or white hole) model for the flow with the equivalent 
event horizon being set down by the condition of the bulk flow 
speed matching the speed of gravity waves ll24l l26l l29l [sTll . 
And so it is that when the conserved spherically symmetric 



accretion flow is decoupled from the space-time geometry in 
the Newtonian limit, one could easily define a metric for an 
acoustic event horizon, which would actually coincide with 
the sonic horizon of the flow. 

This whole aspect of the flow is manifested in its stability 
as well. It has been seen that the stationary inflow solutions 
defined in the Schwarzschild metric are more stable under the 
effects of a linearized perturbation than what they should be in 
the Newtonian construct. Analogously mapped onto the prop- 
erties of the Newtonian limit, this is like the effect of viscous 
dissipation lending greater stability to the flow than what it 
would have been for a perfect fluid. 

Finally, as an interesting aside, it has been shown that the 
steady solution of the relativistic flow leads to the derivation 
of two standard pseudo-Newtonian potentials which are often 
applied to mimic general relativistic effects in the Newtonain 
framework of space and time. 

II. GENERAL RELATIVISTIC EQUATIONS FOR 
SPHERICALLY SYMMETRIC ACCRETION 

A complete description of a conserved spherically sym- 
metric flow in Schwarzschild geometry will require some in- 
dispensible mathematical relations. The first of these is the 
equation for the spherically symmetric line element, which, 
in units of c = 1, can be set down as 

dS^ = ~fdt^ + f-^dr^ +r^dn^ (1) 

with / = f{r,t). This is to be followed by a relation for the 
momentum-energy tensor of a perfect fluid, given by ll32ll 

Tf"" =^ {€+p)v^'v•' +pg^"' (2) 

in which p is the pressure, e is the energy density and is 
the fluid four- velocity, which obeys the relation v^v^^ = — 1. 
In tensorial notation the continuity condition can likewise be 
expressed as ifssll 

ipv"),^ = (3) 

with p being the particle number density of the perfect fluid. 
Some algebra following the use of the two conditions given by 
Eqs. ([T]| and (O will ultimately lead to a modified and explicit 
form of the continuity equation as 

a^(^^^j+^a^(^-) = o. (4) 

This gives one relation connecting the radial flow velocity, v, 
and the local density, p, to each other To solve for each of 
these two fields explicitly, it should be necessary to obtain 
another such relation. This can be derived from the momen- 
tum balance condition. To do so it should be necessary first 
to apply the requirement of energy-momentum conservation, 
T^'•'.^ = 0, on Eq. ©. This wifl give 

(e + p) + v^'v".,) + (e + p) ^ w^w" + g^^y , = 

(5) 
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with e having to be defined separately through the thermody- 
namic relation [33], 
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in which T is the temperature and s is the specific entropy. 

A further definition that is necessary is that of the speed of 
sound, a, which, under isentropic conditions, is expressed in 
terms of the thermodynamic quantities, e and p, as 1.33.] 
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Making use of the foregoing definition, along with the condi- 
tion of constant entropy in Eq. (|6]l, achieved by setting ds = 0, 
and invoking the spherically symmetric line element from 
Eq. ([T]i once again, it should be a straightforward algebraic 
exercise to recast Eq. (|5]l in the form 
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which, incidentally, bears a close resemblance with a similar 
equation derived by Misner & Sharp (32!] in their study of 
gravitational collapse. 

The pressure, p, is connected to the density, p, through a 
polytropic equation of state, p = kp'^, with k and 7 being 
constants, the latter being the polytropic exponent. With the 
help of this equation of state, it becomes easy to show from 
Eqs. (|6]l and (|7J that when the fluid is isentropic, there is a 
relation between a and p that can be written as 
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in which n = (7 — going by the usual definition of the 
polytropic index ||34ll . This connection between a and p now 
makes it possible to see Eq. (O as the second relation, after 
Eq. (HI, that can be written entirely and explicitly in terms 
of V and p. And so with the help of Eqs. @, ^ and (|9]l, it 
should now be possible to establish a complete quantitative 
description of the spherically symmetric flow defined within 
the Schwarzschild metric. 

The stationary solutions of Eqs. (HJi and ([8]l are written as 
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respectively 113311 . In the former solution, the integration con- 
stant, m, is physically the matter flow rate, and p, is the aver- 
age density of particles in the flowing gas. While this solution 
is a direct first integral of the stationary continuity equation, 
the latter solution is not an integral solution. Nevertheless, it 
has certain interesting consequences which become very ap- 
parent in the non-relativistic limit. In this limit, both v'^ <C 1 
and < 1, while r > 2GM. 



Now the static spherically symmetric metric element will 
be given by / = 1 — 2GM/r, from Eq. ([T]i. Therefore, in 
the non-relativistic limit the leading order solution that can be 
obtained from Eq. (fTTT i is 



1 d / 9\ , , dp 

--(.^) + 0'M + -^ = o 



(12) 



in which — f /2f. What Eq. ( fT2] i gives is the stationary 
Euler equation in the non-relativistic limit, with cj) being an ef- 
fective potential driving the stationary flow. At large distances 
(j) behaves like the classical Newtonian potential, but on length 
scales comparable to 2GM, there will be a deviation from the 
Newtonian behaviour. 

Frequently in astrophysics it becomes convenient to dis- 
pense with the full mathematical rigour of general relativity 
and instead, in a Newtonian framework, employ a prescription 
in which general relativistic effects could be represented by 
an effective potential. In this "pseudo-Newtonian" approach, 
many such effective potentials have been suggested according 
to various specific requirements, of which two have been pro- 
posed by Artemova et al. JSi |35|]. Interestingly enough, the 
functional behaviour of these two potentials could be derived 
from the form of (j) implied in Eq. ( fT2] i. Written together, they 



1 - 



2GM 



(13) 



with the latter potential being seen to be actually a special case 
of the former. Both, of course, converge to the Newtonian 
limit on large length scales. 

While all these results can be derived from the stationary 
background flow, to have any understanding of their stability 
under the effect of a linearized time-dependent perturbation, 
it will be necessary to go back to the two dynamic equations 
of the flow given by Eqs. (|4]i and (O. 
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III. LINEARIZED PERTURBATIONS ON STATIONARY 
SOLUTIONS 

For the purpose of carrying out the stability analysis of sta- 
tionary solutions, a standard assumption that is being imposed 
is that the metric is static, i.e. df /dt = 0. This assumption is 
nothing unusual as far as accretion processes are concerned, 
in which the flow is driven by the gravitational field of an ex- 
ternal accretor. This will imply that the gravitational field will 
be unchanging in time 

The perturbation scheme itself will be set down as v{r, t) — 
Vo{r)+v' {r, t) and p{r, t) = pQ{r)+p' {r, t) with the subscript 
"0" indicating stationary values of v and p, and the primes in- 
dicating small time-dependent perturbations about the station- 
ary values. At this point, following the method of Petterson 
et al. fill] and Theuns & David |J2i1 . it will be expedient for 
the perturbative analysis to define a new variable, ij} — pvr^, 
which, as it is very obvious from Eq. (|4]l, is closely associated 
with the matter flow rate, and whose stationary value, ijjQ, as 
it can be seen from Eq. ( fTOl i. is a constant of the motion. The 
first-order fluctuations about this constant stationary value can 
be expressed as 



-0' = {vop' + pov') 



(14) 



Another such relation connecting v', p' and can be derived 



from Eq. (|4|l, and it will read as 



y// + vl dp' pqVq dv' ^ 1 dtp' 



(15) 



With the help of the two foregoing equations it shall now be 
possible to express both p' and v' solely in terms of These 
will be given as 
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respectively. 

The speed of sound, a, is connected to p through Eq. 
and so the perturbation in p has to affect a as well. This is to be 
written as = Oq + (da§/dpo)p'- Once this has been done, 
the linearized first-order fluctuations about the stationary mo- 
mentum balance condition can be extracted from Eq. ^ and 
written as 



f dt dr } po dt dr po ' dr \pQ ' 



(18) 
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Partially differentiating Eq. (fTsT l with respect to time, and spectively, will ultimately deliver a linearized equation of mo- 
making use of Eqs. (fTSI l and ( fTTb to eliminate p' and v', re- tion for tp' as 
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dt 
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dtp' 
dt 
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with the coefficients, ft,"^, having to be read from 



(f + vl-vlal) 



At this stage it would be very much instructive to exam- 
ine the features of Eq. (fT9] l in the non-relativistic limit, where 
both Vq and Oq are vanishingly small compared to unity while 
/ itself assumes the value of unity. As a consequence, in the 
non-relativistic limit it will become eminently possible to re- 
duce Eq. iT% to a compact form given by 



do. {h^^dptP') = 



(20) 



and 



with the Greek indices a and /3 running over t and r, as they 
did for the fully relativistic case earlier 

In Lorentzian geometry the d' Alembertian for a scalar field 
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in curved space is given in terms of the metric g^^, by ifToll 

A^= -La^v^g^'-a,^) (21) 

with being the inverse of the matrix impHed by g^^. Com- 
paring Eq. ( l20l ) with Eq. (l2Tl i. it would be tempting to look 
for an exact equivalence between h"^ (in the non-relativistic 
limit) and \/— gg^"- This, however, cannot be done in a gen- 
eral sense. What can be appreciated, nevertheless, is that 
Eq. ( I20I 1 gives a relation for ijj' which is of the type given by 
Eq. (I2TI ). The metrical part of Eq. ( |20] |. as given by the values 
of ft,"'^, may then be extracted, and its inverse will incorpo- 
rate the notion of a sonic horizon of an acoustic black hole 
when Vq — a^. This point of view has some features similar 
to the metric of a wave equation for a scalar field in curved 
space-time, obtained through a somewhat different approach, 
in which the velocity of an irrotational, inviscid and barotropic 
fluid flow is first represented as the gradient of a scalar func- 
tion, and then a perturbation is imposed on this scalar func- 
tion Lia 23J. 

While all this similarity is undoubtedly pleasing to note in 
the non-relativistic limit, it is also quite obvious from Eq. ( fT9] l 
that in the fully general relativistic treatment there is a break- 
down of the symmetry that leads to the devising of an ana- 
logue gravity model for the fluid flow. This is at striking vari- 
ance with the conclusions arrived at by Bilic ll20ll . It is not 
very difficult to discern that the difference arises because of 
the way the perturbative studies have been prescribed in the 
two cases — the scalar potential approach of Bilic ll20ll . and 
the continuity equation approach in the present case. 

Having said so, it must nevertheless be pointed out that in 
some respects at least the fluid analogue approach is not en- 
tirely lost even for the general relativistic flow being studied 
here. Speaking by analogy, in the shallow layer flow of a 
perfect liquid (without any viscosity) it has been shown that 
an analogue black hole or white hole model is very much 
a fact |2l|] . Viscous dissipation in the flow, on the other 
hand, adversely affects the invariance that makes the ana- 
logue model possible, even as it helps the hydraulic jump phe- 
nomenon itself to happen ll24l[3lll . And in a like fashion while 
it is easy to recognize the concept of fluid analogue gravity 
for spherically symmetric accretion in the non-relativistic sit- 
uation, the coupling of the geometry of space-time with the 
perturbed field in the fully general relativistic scenario acts in 
the manner of a dissipative effect. At least through the con- 
tinuity equation approach this precludes any hope of building 
an acoustic black hole model within the general relativistic 
flow, but at the same time this has a more favourable role to 



play as regards the stability of the stationary solutions. 



IV. STABILITY ANALYSIS : STANDING WAVES 

In treating the perturbation as a standing wave it will be 
essential to identify proper boundary conditions to constrain 
the wave at two spatial points. Between these two points the 
properties of the standing wave could be studied. One bound- 
ary condition for the perturbation is conveniently fixed at the 
outer boundary of the stationary flow itself, where any solu- 
tion naturally decays out. So should any perturbation imposed 
on it. 

In determining an appropriate inner boundary condition, on 
the other hand, one encounters much greater difficulties. The 
nature of the accretor itself has an influential role to play in 
this matter If it is a black hole then the infalling matter has 
to cross the event horizon at the maximum possible rate ll36ll . 
and the only feasible solution should be a transonic one, which 
will smoothly pass through a singular point in the flow |37*]. 
On the other hand, if the accretor is a compact object with a 
physical surface, then the inner boundary condition becomes 
subject to many complications, which in turn will leave its im- 
print on the character of an inflow solution in the vicinity of 
the stellar boundary. The solution could either be transonic or 
shocked (discontinuous) subsonic or continuously subsonic. 
These distinct aspects have been discussed at length by Petter- 
son et al. [11] and Theuns & David |,12.] . In this context Mon- 
crief 0] has pointed out that as long as the stellar surface is 
enclosed within the sonic surface, his perturbative arguments 
hold both for a black hole as well as a compact star. 

The standing wave analysis that has been pursued here will, 
of necessity, require the background solution to be continuous 
everywhere and globally well behaved. Besides this the wave 
will also have to die out at the two chosen boundaries. Now 
the only solutions that will meet all these requirements in a 
general sense are the purely subsonic solutions. While these 
may not entirely be representative of the precise manner of 
infall in the general relativistic scenario, a mathematical study 
of their stability will reveal the true extent of the influence 
that the coupling of the flow with the geometry of space-time 
will have vis-a-vis what it is in the Newtonian flat space-time 
limit. Stability of subsonic flows in the latter situation has 
been studied thoroughly by Petterson et al. ifllll and Theuns 
& David fH . 

The solution ijy{r, t) = p{r) exp {—\u)t) is to be applied in 
Eq. (fT9T l and the resulting expression multiplied throughout by 
p. This will lead to 



H'V.' + i { A (,,.y ) _ ,.viL [m + + [, (/ + ,5)] - ± . ,22) 

I 



which will then have to be integrated over the entire spatial 



range within which the standing wave is continuously dis- 
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tributed. At the two boundaries of this range the amplitude 
of the wave is required to vanish. So the "surface" terms 
obtained from integrating Eq. (l22l i will also have to vanish at 
the boundaries. This will leave a residual quadratic equation 
which will be a dispersion relation for w. This relation will be 
in the form 



Au"^ - 2iBlo + C = 
with the three coefficients above being read as 



(23) 



is otherwise a conservative system. For a conserved flow in 
the Newtonian regime, Petterson et al. have shown that 
the perturbation will have a constant amplitude. Any decay of 
the amplitude of the standing wave could only be reproduced 
when one accounts for viscosity in the flow [13,1 . Coming 
back to the conserved general relativistic case, the only pos- 
sible explanation for the decaying behaviour can be that the 
coupling of the flow with the geometry of space-time acts in 
the manner of an "effective" dissipation. And consistent with 
this line of thinking it can also be shown that in the Newto- 
nian limit one does indeed regain the expected constancy of 
the amplitude of the standing waves. 



B = 
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[ln(/ + t;2)]dr 



and 



C 



dp d 



f + Vq dr dr 



[P{f + vl)]dr, 



respectively. 

Under the condition that for inflow solutions i; < and 
(dpo /dr) < 0, it shall be easy to verify that {B /A) < by 
referring to Eq. ( fTTI ). Therefore, the discriminant of the so- 
lution of Eq. ( l23l l will hold the key regarding the stability of 
the standing wave. Once again, for subsonic solutions, it can 
be argued that {C/A) < 0. And so, if \C/A\ > {B/Af, it 
will imply that the time-dependent part of the standing wave, 
given by exp(iLjt), will have an oscillatory nature with the 
amplitude of the oscillation being damped in time. Hence 
the background solution will be stable. If, on the other hand, 
\C/A\ < (B/A)^, then there will be two real roots of w, both 
negative, indicating that the perturbation will be overdamped. 
So, one way or the other, the amplitude of the standing wave 
will be damped, lending stability to the stationary background 
solution in strong measure. 

This is a rather intriguing state of affairs indeed. The de- 
cay of the amplitude of the standing wave would imply that 
something in the nature of a dissipative effect is active in what 



V. STABILITY ANALYSIS : TRAVELLING WAVES 

The manner in which the stability of the flow is influenced 
by its coupling with the space-time metric could also be ex- 
amined by fashioning the perturbation to be a high-frequency 
travelling wave. A comparison could then be made with the 
corresponding analysis carried out in the Newtonian structure 
of space and time by Petterson et al. [11] who argued that the 
travelling waves could cause a growth in the fluctuations on 
the flow rate, but would not drive the background flow from 
it stationary profile. To this extent the stability of the flow 
should be preserved. 

High-frequency travelling waves are to be first defined pre- 
cisely in the present context by the fact that their wavelength 
should be much smaller than the Schwarzschild radius of the 
black hole. This will imply that the frequency, uj, should be 
correspondingly large. With this restriction on uj, the spatial 
part of the perturbation can then be prescribed in terms of a 
power series as 



Pu;{r) = exp 



J2 ^"'^'(^ 



.i=-i 



(24) 



and this is then to be applied to a slightly modified rendering 
of Eq. (I22I 1 which goes as 



dr^ 



dh" 
dr 



dp 
dr 



dr dr J 



From the result of this extended algebraic exercise all coef- 
ficients of cj^ are to be collected and their sum is to be set to 
zero. This will give a solution for /c_i, which will look like 



k-i 



(h^y 



dr. (26) 



Similarly summing up all the coefficients of uj to be zero, and 
applying the value of as Eq. d26] l gives it, will deliver a 



solution for fcn as 



\J{h*''f - h"ht^ 



1/2 



(27) 



Likewise, the solution of fci could be found by setting the sum 
of the coefficients of to be zero. In terms of fc_i and fco 
this can be expressed as 
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dr dr 



\n{f + vl)]=Q 



(28) 



For reasons of self-consistency it shall be necessary at this 
stage to show that successive terms in the power series given 
by Eq. (l24l l will obey the requirement that a;^'|/c/(r)| ^ 
|fc;_i_i(r)|, which will also imply that the power se- 
ries will converge rapidly with increasing I, and so it can be 
truncated after the first few terms. Mindful of the fact that lo is 
large, the first three terms, involving fc_i, fco and ki, respec- 
tively, conform to this self-consistency requirement. A simple 



J 



f + v'o 



asymptotic check suffices to show fc_i ^ r, fco ^ Inr and 
fci ~ r~^. In any case both fc_i and ki make contributions to 
the phase of the travelling wave. In this linearized treatment, 
therefore, the most conspicuous contribution to the amplitude 
comes from fco, and any impression of the stability of the flow 
can be unambiguously derived from this term only. 

Using only the solutions of fc_i and fco, the dominant prop- 
erties of the perturbation could be set down as 



1/2 ^ 



dr > e 



(29) 



with ^ being a constant, and with the positive and negative 
signs, placed together, indicating a superposition of incom- 
ing (corresponding to the negative sign) and outgoing (corre- 
sponding to the positive sign) travelling waves, respectively. 

It will now be very much worthwhile to scrutinize the ex- 
pression for fco more closely and see how the geometry of 
space-time makes its contribution to stability. Using the de- 
rived values of /i**, /i*'' and /i^*" it can be shown that fco picks 
up a term that goes as In /. By virtue of the fact that / < 1, 
the logarithm of / will be negative, and consequently its ef- 
fect on fco would be to detract from its value in the Newtonian 
limit (where / — 1). And so where general relativistic effects 
will have to be accounted for, the amplitude of the travelling 
waves will become subdued. Once again it will not be difficult 
to see that this effect is owed entirely to the coupling between 
the flow and the curvature of space-time. 

To delve into some more details, the amphtude of the per- 
turbation could be recast very simply as 



1/4 



(30) 



a form which is quite helpful in shedding a clear Ught on the 
asymptotic behaviour. In the Newtonian limit, it is obvious 
that I?/;'! ~ (ao^o)"^^^, a state of affairs whose stability has 
been cogently argued for by Petterson et al. An equal 

measure of stability is to be seen near the event horizon as 
well, with a smooth passage for the travelling wave through 
the sonic region somewhere in between. After this it becomes 
easy to argue for the asymptotic stability of the steady solu- 
tions. 



VI. CONCLUDING REMARKS 



Some parting comments would well be in order It has been 
shown already that the manner in which the perturbative study 
has proceeded, beginning with the continuity condition of the 
flow, has led to a failure in establishing an acoustic geometry 
for general relativistic spherically symmetric accretion. On 
the other hand, as far as stability is concerned, this approach 
has been in perfect qualitative conformity with earlier studies 
which have dwelt on the question of the stability of the flow 
solutions. Indeed as regards stability in particular, it has been 
argued and shown that general relativistic effects enhance the 
stability of the stationary solutions. This fact actually opens 
up an interesting possibility. 



Schwarzschild space-time defines a geometry of closed cur- 
vature. A standing wave in this geometry exhibits a damp- 
ing of its amplitude, and so a stable behaviour is implied. In 
the Newtonian limit, the standing waves continue to have a 
constant amplitude. In this respect the behaviour may once 
again be viewed to be stable. However, the situation could 
become radically different in the geometry of open curvature. 
One might conjecture that in this case, simply because of the 
nature of the geometry, there will be an unstable behaviour, 
manifested through a growth in the amplitude of a standing 
wave. Indeed, this speculation is not without its foundations. 
In a different context, but with good consonance, similar fea- 
tures are to be seen in CMB anisotropy in compact hyperbohc 
spaces iMm. 
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